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O ABSTRACT 

A factorization theorem for P-wave quarkonium production, recently derived by Bodwin, Braaten, 
Yuan and Lepage, is applied to T — > XcJ + X, where XcJ labels the 3 Pj charmonium states. The 
widths for XcJ production through color-singlet P-wave and color-octet 5*- wave CC subprocesses are 
computed each to leading order in <y, s . Experimental data on T — > J / ifj + X is used to obtain an 
\Q . upper bound on a nonperturbative parameter (related to the probability for color-octet o -wave CC 

hadronization into P-wave charmonium) that enters into the factorization theorem. 



Factorization theorems play a basic role in perturbative QCD calculations of 
many hadronic processes. A well known factorization theorem for the decay and 
production of S-wave quarkonium follows from a nonrelativistic description of heavy 
quark-antiquark (QQ) binding [I]. Nonperturbative effects are factored into Rs(0), 
the nonrelativistic wave function at the origin, leaving a hard QQ subprocess matrix 
element that can be calculated in perturbation theory. 

Remarkably, the correct factorization theorems for the decay |2j and produc- 
tion H of P-wave quarkonium have only recently been derived. These new theorems 
resolve a long standing problem regarding infrared divergences which appear in some 
cases to leading order in the rates for P-wave QQ states Q. In previous phenomeno- 
logical calculations, the divergence was replaced by a logarithm of a soft binding scale 
[H, However, a rigorous calculation requires that one consider additional compo- 
nents of the Fock space for P-wave quarkonium, such as \QQg), where the QQ pair 
is in a color-octet S-wave state, and g is a soft gluon 0, |3|. 

A renewed study of the decay and production of P-wave quarkonium is there- 
fore of considerable interest, since one may gain new information on a nonperturba- 
tive sector of QCD that has largely been neglected in the quark model description 
of heavy quarkonium. This is also of practical consequence; for example, J/i/j pro- 
duction provides a clean experimental signature for many important processes, and 
P-wave charmonium states have appreciable branching fractions to J/i/j. 

* Presented at the Workshop on Physics at Current Accelerators and the Supercollider at Argonne 
National Laboratory, June 1993. 



In this paper the factorization theorem for P-wave quarkonium production is 
applied to T — ► Xcj + X, where Xcj labels the 3 Pj charmonium states. The widths for 
Xcj production through color-singlet P-wave and color-octet ^-wave cc subprocesses 
are computed each to leading order in a s . Experimental data on T — > J/ijj+X is used 
to obtain an upper bound on a nonperturbative parameter (related to the probability 
for color-octet S"-wave cc hadronization into P-wave charmonium) that enters into 
the factorization theorem. The bound obtained here adds to the limited information 
so far available on the color-octet mechanism for P-wave quarkonium production. 
The color-octet component in P-wave decay was estimated in Ref. || from measured 
decay rates of the xa an d Xci- A rough estimate of the color-octet component in P- 
wave charmonium production was obtained in Ref. || from data on B meson decays; 
however, an accurate determination in that case requires a calculation of next-to- 
leading order QCD corrections to the color-singlet component of B — > Xcj + X, which 
is so far unavailable ||. 

The factorization theorem for P-wave quarkonium production has two terms, 
and in the case of T decay takes the form: 

T(T^Xc.j + X) = Hifi(T^cc( 3 Pj)+X;ri 

+ (2 J + l)H'Mt 8 (T -> cc( 3 S 1 ) + X). (1) 

f i and f § are hard subprocess rates for the production of a cc pair in color-singlet 
P-wave and color-octet S'-wave states respectively. The quarks are taken to have van- 
ishing relative momentum. The nonperturbative parameters Hi and H 8 are propor- 
tional to the probabilities for these cc configurations to hadronize into a color-singlet 
P-wave bound state. Hi, H 8 and Tg are independent of J. This factorization theorem 
is valid to all orders in a s and to leading order in v 2 , where v is the typical center- 
of-mass velocity of the heavy quarks. The hard subprocess rates are free of infrared 
divergences. r\ and H 8 depend on an arbitrary factorization scale [i in such a way 
that the physical decay rate is independent of fi. In order to avoid large logarithms 
of mr//-f in T 1 , \i of O(ra-r) should be used. 

Hi can be expressed in terms of the P-wave color-singlet cc wave function, 
Hi pa 9|Pp(0)| 2 /(27rm^) ph 15 MeV, where the numerical estimate was obtained in 
Ref. @ from measured decay rates of the Xci and Xc2- H' 8 cannot be rigorously 
expressed perturbatively in terms of Rp, since it accounts for radiation of a soft 
gluon by a color-octet cc pair. The scale dependence of H' 8 (/j,) is determined by a 
renormalization-group equation 0, || which (at leading order in a s (/•*)) gives 0: 



H' 8 (m b ) = H'M + 




Hi (2) 



(for /i < m c ), where (3 n = (33 — 2n)/6. If H' 8 (n ) is neglected in the limit of large m& 
one obtains H' 8 (mb) ~ 3 MeV, using a s (/i ) ~ 1 §. While one might not expect the 
physical value of m b to be large enough to neglect H 8 (fi ), an estimate for H' 8 {m b ) 
obtained in Ref. J| from experimental data on B meson decays is consistent with 
the above result. 



A calculation of f i and fg in Eq. ([[[) each to leading order in a s can be 
obtained from a calculation of the infrared divergent width for T — ► cc( 3 Pj)+ggg, 
where the cc pair is in a color-singlet P-wave state: 



r d iv(T -> cc^Pj) + ggg; /x ) 
20a 5 s Gj 



3 7 7r 3 m x L 



^u(//) + (2J + l 




Pi- (3) 



T\j and JF 8 are dimensionless infrared-finite form factors, /xq is an infrared cutoff on 
the energy of soft gluons, and /x is an arbitrary factorization scale [the /x dependence 
of T\j exactly cancels that of the explicit logarithm in Eq. @]. The constants in 
Eq. (H) include a color-factor of 5/216 and phase space factors, including 1/3 for T 
spin-averaging, and 1/3! for the phase space of the three indistinguishable gluons [cf. 
Eq. (§) below]. Gj is related to the usual S"-wave bb nonrelativistic wave function, 
Gj ~ 3|-Rj(0)| 2 /(27rm^) pa 108 MeV, where the numerical value is obtained from the 
electronic decay rate of the T [|]. 

The hard subprocess rates of Eq. ([!]) are identified from r^v by using the 
perturbative expression for the infrared divergence in Pg, obtained by neglecting the 
running of the coupling ||: Pg(/x) ~ (16/277r)a s ln(/x//Xo)Pi. Thus: 

2Drv 5 G T 

r\(T - ce( 8 Pj) + ggg; /x) = ^j-^M, (4) 



and 



9Dry 4 G T 

f 8 (T - ccfSi) + gg) = ^ri^s- (5) 

6' it 6 m x 



Note that f i is suppressed by 0(a s ) compared to f 8 . However, the nonperturbative 
parameters Pi and P 8 which accompany these subprocess rates in Eq. (JJ) are inde- 
pendent, hence a s Hi need not be small compared to P 8 0, |J. We therefore proceed 
to calculate r\ and r 8 each to leading order; all further corrections to P-wave pro- 
duction are then guaranteed to be suppressed by at least one power of a s compared 
to what is included here. 

In order to extract T\j and T% individually, it is necessary to explicitly identify 
the infrared logarithm in the calculation of Y^ y . This can be done analytically, as 
described in the following. There are 36 0(a 5 s ) diagrams contributing to Tdiv One of 
these is shown in Fig. |l|. Define the invariant amplitude Aij(2, 3; 1) corresponding to 
the sum of all Feynman diagrams where gluon "1" is radiated from the charm quark 
line. The amplitude is readily computed using expressions for S- and P-wave QQ 
currents given in Ref. 

m (2 o. 11 = m T m x P,(2, 3)^(1) 

A ' ' ' l(k 2 + h) ■ k 3 ] i(k 3 + k 4 ) ■ k 2 ] [(k 2 + h) ■ k A ] k\ (k ■ hy ' 1 1 



^Overall factors in the quark currents including couplings, color amplitudes, and wave functions have 
been accounted for in Eq. (H). 



Figure 1: One of the 36 0{a\) diagrams contributing to T — > Cc( 3 Pj) + ggg. 



where 

e%B^(2, 3) = {e 4 • e 2 [-£;4 • fae 3 • fato ■ fa - fa ■ k 3 e 3 ■ fae -fa- fa- k 3 k 2 ■ k 3 e ■ e 3 ] 
+ eo • e 3 [fa ■ fa^ • fati -fa + fa- k 3 t 2 ■ fae 4 -fa-fa- k 2 e 4 ■ k 3 e 2 ■ k 3 }} 
+ {2 <-> 3} + {3 <-> 4} (7) 

(eo is the polarization of the T), and 

^C^ =0 (l) = \J\[ti- e 4 fa ■ fa - ei • /c 4 e 4 • fa] (m 2 x + fe • fa - kfj , 

e 4/U C^ =1 (l) = \m x k\e a ^e a e{e[k[, (8) 

e 4M C^ =2 (l) = sj\m\ (fa ■ fae°tei + k%k{e x ■ e 4 - fae{e x ■ fa - fc?efe 4 ■ fa) e af3 . 

e a is a spin-1 polarization vector and e a/3 is a spin-2 polarization tensor. For conve- 
nience the virtual gluon is labeled in Eqs. by polarization e 4 and momentum 
fc 4 (fa = P — fa — k 3 = k + fa). Terms which vanish due to the on-shell conditions 
€i ■ fa = (i — 1, 2, 3) and eo • P = have been dropped. 
The overall factors in Eq. ([3]) are such that: 

*,00 + <2J + 1)£ - 

3/4*i] £ [M 2 ,(2,3;1) + 2A1,,(2.3;1)A<.,(1,3;2)], (9) 

spins 

where $ n denotes (infrared-cutoff) n-body phase space, normalized according to 

$ n [p - P1 , . . . , Pn ] = j n ^5 4 (p - ( io ) 



The factor of 3 on the right hand side of Eq. (|9|) accounts for symmetrization of 
Aij(2,3; 1) under gluon label interchanges 1 2 and 1 •*->• 3, taking account of the 
symmetry in the three gluon phase space. 

The infrared divergence comes entirely from the first term in square brackets 
in Eq. (^), and is due to the P-wave charm quark propagator l/(k ■ ki) 2 in Eq. (Eft. 
It is therefore advantageous to organize the four-body phase space integral in Eq. (pi) 
by taking the invariant mass of the Xcj and gluon "1" as one integration variable || 

d[* 4 ] = / d(k 2 23 ) / d{k\ x ) 

X $ 2 [P -> k 23 , k lx ] $2^23 -> &2, *&] -> fcl, fc], (11) 

where m| 3 = k 23 . Note the infrared cutoff \xq on the energy of gluon "1" in the rest 
frame of the \ c j. 

The infrared logarithm on the right-hand side of Eq. @ can now be identified 
analytically by observing that P M (2, 3)Cj(l) in Eq. (|j) is given by a sum of terms 
each containing exactly one factor of k±, if k± = P — k% — kg is used to eliminate the 
virtual gluon momentum. With this convention, one has 

E ^2,3;l) = ^^M , (12) 



spins 



(k ■ kij 



where k\ appears explicitly in the function ^j{k\\ P, k, k%, k 3 ) only in the combination 
k\jk ■ k\. JF 8 is then given in terms of a manifestly infrared-finite three-body phase 
space integral, taking account of the fact that $2(&ix — ► ki, k) = ^k ■ k\/k\ x j dQ{ x , 
where f2* x is the center-of-mass solid angle of the two body system: 



Z/7T f \ rri T ~ m x) 

(2J+m = 32^- y o d(ki 3 )^ 2 [p^k 23 ,k] 

x$ 2 [k 23 ^k 2 ,k 3 } I dn* lx 7j(huP,k,k 2 ,k 3 ), (13) 

where 

h 

fci = lim — . (14) 
fc-fci^o k-ki 

The finite four- vector k\ is readily expressed directly in terms of k 23 and . An 
expression for T\j can be obtained from Eqs. ([]) and ( |T3| ) by analogy with the 
identity f dxf(x)/x = f{0)kix + J dx\f(x) - f(0)]/x. 

The contraction of currents and sum over polarizations in Eqs. (||) and (|9]) 
were performed symbolically using REDUCE [0] (leading to lengthy expressions, par- 
ticularly for J = 2). The Xcj spin sums were done using (see e.g. Ref. M): 

k k 

E__ , rh fl nj L' _ q~) 

m\ 

e ^e a p = \ \P m V v j3 + V^V ua ] - lVfiuVa/3- (15) 



Figure 2: Color-octet form factor J-'g as a function of m x /mx- 



The phase space integrals were evaluated numerically using VEGAS [ I0| ; modest 
integration grids are found to give very good convergence. The fact that T% should 
be independent of J provides a stringent check of these calculations, given that the 
three currents Cj have very different structures [cf. Eq. (0)]. This was verified 
explicitly in numerical calculations of Eq. (|i"3l), to better than a few tenths of a 
percent for all m x /m-x on a modest integration grid. Figure shows the numerical 
results for JF§ over a range of hypothetical meson masses. In Fig. ||| results for Tij{ix) 
are shown using a factorization scale \x = rax- 

The available experimental data on charmonium production in T decay is for 

the J 

= (1.1 ± 0.4) x 1(T 3 CLEO [11], 
B cxp (T -> J/ip + X) \ < 1.7 x 1(T 3 Crystal Ball 0, (16) 

< 0.68 x 10~ 3 ARGUS @. 



An upper bound on H' % can be extracted from this data by computing the "indirect" 
production of J /if) due to the Xcj states. Assuming that radiative cascades from Xci 
and Xc2 dominate, with branching fractions -B ex p(Xci — ► lJ/4 1 ) ~ 27% and -B exp (Xc2 — >■ 
7 J /if)) ~ 13% ||, the results presented here give: 

8V 7 \ 2.9 x 10- 5 J V 7 

The first number in brackets above comes from the color-octet subprocess rate Tg, 
and the second number from the color-singlet rate T\. The experimental value for 
the total width r to t(T) w 52 keV was used, along with a s (m r ) w 0.179 QTJ] , and 
the values of i^i and given above. 



Figure 3: Color-singlet form factors T\j as functions of m x /fftx : J = (short-dashed line), 
J = 1 (long-dashed line), J = 2 (solid line). The form- factors were evaluated using a factorization 
scale n = my. 



Equation ( |TTD yields the bound ifg(m T )<25 MeV using the ARGUS upper 
limit, which is consistent with the other measurements. This bound is considerably 
larger than an estimate H' 8 (m b ) « 3 MeV based on B meson decays although a 
calculation of next-to-leading order QCD corrections to the color-singlet component 
of B — * Xcj + X is required before an accurate determination of H' % can be made in 
that case ||. 

This raises the possibility of significant direct production of J/ip in the decay 
of the T, unless the branching fraction turns out to be considerably smaller than the 
ARGUS bound. Mechanisms for direct T — > J/ip + X in perturbative QCD were 
first discussed in Refs. fl4| and Hl5| . The direct production rate is suppressed by 



0(a 2 ) compared to the P-wave color-octet production mechanism considered here. 
However, the nonperturbative matrix elements which enter into P-wave production 
are of 0(v 2 ) relative to the corresponding parameter for S'-wave production. 

The full 0{a & s ) perturbative QCD amplitude for direct T — ► J ftp + X was 
recently evaluated in Ref. |16[ , corresponding to one-loop diagrams for T — > J/ijj+gg, 
and tree diagrams for T — > J/ip + gggg. The 0(ct 2 s a 2 ) electromagnetic amplitude for 
the two gluon decay mode was also evaluated. Unfortunately, only a crude estimate 
of the required phase space integrations was made in Ref. |16| (there is a costly 
convolution with a numerical calculation of the loop integrals for T — > J/ip + gg). 
Nevertheless, the calculation of Ref. |R| suggests a branching fraction for direct 
production of a few x 10~ 4 . This would lead to a considerable reduction in the bound 
on ifg extracted from Eqs. ( Jl6| ) and ([I] 



^H' s (fi) increases by only f» 0.3 MeV in the evolution from [i = rrib to fl = mr [cf. Eq. (0)]. 



To summarize, a complete calculation was made of the leading order rates 
for T — > Xcj + X, through both color-singlet P-wave and color-octet S'-wave cc 
subprocesses. Experimental data on J/if> production was used to obtain an upper 
bound on the nonperturbative parameter H$, related to the probability for color-octet 
S'-wave cc hadronization into P-wave charmonium. Improved experimental data, and 
a definitive calculation of the direct J /if) production rate along the lines of Ref. fl6|l , 
would allow for an accurate determination of H' % from the results presented here. 

I am indebted to Eric Braaten for suggesting this problem, and for many 
enlightening conversations. I also thank Mike Doncheski, John Ng, and Blake Irwin 
for helpful discussions. This work was supported in part by the Natural Sciences and 
Engineering Research Council of Canada. 
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